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SUMMARY 

This paper presents some recently added capabilities to RAGGS, Rockwell Automated Grid Generation 
System. Included are the trimmed surface handling and display capability and structured and unstructured 
grid generation on trimmed Bezier and NURBS quilted surfaces. Samples are given to demonstrate the new 
capabilities. 


INTRODUCTION 

As described in Ref. 1, RAGGS (Rockwell Automated Grid Generation System) provides an integrated 
CFD (Computational Fluid Dynamics) environment where a user can 1) import a CAD (Computer Aided 
Design) generated configuration or a geometry from other sources, 2) build a CFD grid directly on the im- 
ported geometry, 3) set up flow solver input, 4) postprocess the flow solution, and 5) display the flow solution 
data and geometry simultaneously. The program structure of RAGGS is depicted in Fig. 1. It consists of a 
series of modules to handle various tasks in the CFD process: starting with a given geometry, usually input 
as an IGES (Initial Graphics Exchange Specification) file, through construction of CFD geometry database, 
domain blocking, surface grid generation, volume grid generation, flow simulation, to postprocessing of the 
flow solution to obtain data useful for engineering design and analysis. A file containing the geometry data 
read by RAGGS is referred to as a RASCAL (Rockwell Automated System for Computer Aided Lofting) file. 
As shown in Fig. 1, each module is connected to the common database shown in the central hub. Through 
the GUI (Graphical User Interface) module and/or grid file module, RAGGS enables user to quickly com- 
municate with some state-of-the-art flow solvers, flow analysis software, and grid generators as listed in the 
figure. The RAGGS surface grid generation is operated directly on the quilts of rational Bezier patches 
and/or NURBS (Non Uniform Rational B-Spline), each of arbitrary order. This minimizes the interaction of 
CAD/CFD groups and eliminates the need of re-constitution of the imported geometry, which often requires 
hundreds of hours of manipulation to achieve a data format suitable for the grid generation software package. 

This paper presents some recently added capabilities to RAGGS. These include trimmed surface/ curve 
handling and display, unstructured surface grid generation, and structured/unstructured grid generation 
on trimmed Bezier and NURBS quilted surfaces. The materials presented in this paper are intended to 
complement that given in Ref. 1. 


GEOMETRY DATABASE 

The coordinates of a surface in space can be described in terms of two parametric variables. A patch is 
a mapping from a space defined by these two parametric variables to the three Cartesian coordinates in 3-D 
space. The untrimmed domain of each patch is the rectangular region in parameter space, so < s < Si and 
to < t <ti, where usually s 0 = t 0 = 0 and s x = ti = 1. Additionally, by defining active and inactive regions, 
trim curves may be used to further restrict the domain of a patch. These trim curves are nonintersecting 
closed loops of Bezier curves or NURBS in the parametric space of a patch. 

For a general configuration, aircraft or otherwise, there is no uniform parameterization which can describe 
the entire surface. However, the surface can be described by a set of individually parameterized patches, 
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possibly with trim curves. Any subset of these patches, having surface continuity, can be said to constitute 
a quilt. Any configuration or any number of configurations can be completely described by a number of such 
quilts. RAGGS uses such quilts of rational Bezier patches and/or NURBS, each of arbitrary order, to create 
grids which lie exactly on the parameterized surface. These patches may have 3 or 4 sides and need not be 
distributed in any type of uniform manner. 

Rational Bezier Patch and NURBS 

Bezier and NURBS curves in 3-D space are composed of a point which varies parametrically with a single 
parameter. Surfaces are composed of a curve which itself varies parametrically in 3-D space according to 
another parameter. Consider a such NURBS curve, of degree n, parameterized in t. The curve may be 
written in terms of the B-Spline basis functions as follows: 





vj — the weight coefficients 

gj — the control points 

( &"(<) — the B-Spline basis functions 


( 1 ) 
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A knot vector is defined by the nondecreasing sequence t_„, . . . , to, U, • • • > ^K n > t K n + 1 • The B-Spline basis 
functions are defined by the following recursion, 
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The basis functions have the property that 
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For a rational Bezier curve, K n = n, is equal to the degree of the curve, the parameter domain becomes 
[0, 1], the knot vector becomes J_ n +i , . . . = 0, an< i • • • > *n>^n+i — an< i basis functions become 
the Bernstein functions, 

(4) 

If we let the weights, vj, and the control points, g jy be interpolated functions of another parameter s, 

Km Y.w Jx p Jt b^(s) 

Vj(s ) — ^2lVjib™(s) Qi(s) ~ (5) 


»— 0 


Yjh jib?(s) 

i=0 

Then Eq. (1) together with Eq. (5) constitutes a surface patch parameterized in s and t as follows: 
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with pij being the control points and Wij the weights. 


For fixed i, Eq. (6) may be written as 


Km 

£ 

i=0 


H(s) = P(s,t) = 


Kn 


Kn E (0 



J=0 E «iW) 

i=o 




i=0 


/Cm 

E 

i=o 


E < m ?(0 

J=0 


Kn 


K(») 


E«.-(0 fc r(«) 

t=0 


with 


/C« 


Kn 

«.(*) = E "#*?(*) 

1=0 


E“i<ft« 6 ?(0 


^■(0 = 


j=0 


Equations (7)-(8) define another set of curves parameterized in s for fixed t, i.e., for fixed t, 
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with U{ being the weights, h{(t) the control points, and the B-spline basis functions. 


In RAGGS, all the surface grid generators are formulated in the 3-D physical space, i.e., the grid points 
are computed directly in terms of the three Cartesian coordinates, x, y, and z. These space points are 
then placed exactly on the nearest quilted surface point by iterating for the patch number and parametric 
variables. A key element in a formulation using discrete patches is the boundary data, which is a description 
of how the parametric variables change as one crosses from the boundary of one patch to the boundary of 
an adjacent patch. 


Definition Of Active Region Of Patch 


As mentioned before, a surface is descibed by a quilt of Bezier patches or NURBS, possibly with trim 
curves. Each patch is parameterized in s and t. The untrimmed domain of each patch is sq < s < s\ and 
to < t < tu where usually s 0 = f 0 = 0 and = t\ = 1. Trim curves (if any) are nonintersecting Bezier 
curves or NURBS forming closed loops in the parametric space of a patch. 



i = 


( 10 ) 


Here, Ci(or) is the ith parameterized trim curve, a is the parameter, and n is the number of trim curves. 


The entire domain of an untrimmed patch is defined to be active. Patches which have trim curves have 
both active and inactive regions. By convention, the active region of a trimmed patch is the region that lies 
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on the left of each trim curve as a point advances along the curve. Here, left of a trim curve at a point C t (<r) 

is the direction of the vector formed as the cross product (in the order specified) of the surface normal and 

the tangent vector to the trim curve at Ci(c r). 

Patch Boundary Data 

The patch boundary data is a description of how the parametric variables change as one crosses from 
the boundary of one patch to the boundary of an adjacent patch. Boundary data for each patch is given in 
terms of a series of parameter ranges for each side. Figure 2 shows a typical boundary connectivity between 
two adjacent patches (in this case, patches p and 9). For patch p, there may be data for ranges, 

S bl < S < Sfcj 

^63 — & ^ ^64 (H) 

For each parameter range, an adjacent patch number and a connectivity matrix of six numbers are given. In 
Fig. 2, we assume that a point has crossed the boundary at (s&,f&) with 

h = t 1, s 0 <s b <si (12) 

Then the boundary data for the side t b = h and the range which includes s b will give the new patch number 
q and the corresponding parameter values of the boundary point (t b ,s b ) on this patch q as follows: 

('* )*,**, 

where Ojj are the transformation matrix coefficients. 

Trim Curve Boundary Data 

Consider the trim curves of Eq. (10). Analogous to the patch boundaries, boundary data for trim curves 
is given for a series of parameter ranges of the curve being crossed as follows: 

<76, < o < cr&j 

a b3 < a < a bi (14) 
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patch p 


When crossing trim curve boundaries, the boundary data which includes the parameter range of the crossing 
point < 7 b will give the new patch and the corresponding trim curve on the new patch, and the parameter 
value of the boundary point on the corresponding trim curve as folows: 
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where ai and a 2 are the transformation matrix coefficients. 
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Geometry Creation and IGES/RAGGS Interface 

The basic method of geometry input is through IGES files created by any CAD system. The data in 
the IGES file must be translated into surfaces in the form of NURBS and Bezier Patches. These are the 
fundamental surface forms used by RAGGS. Translating the data into the RAGGS internal format may take 
a number of different paths, depending on the type of data in the IGES file. A flow diagram of the process is 
shown in Fig. 3. The processing of the IGES file yields four different types of information which are stored 
for later use. 

NURBS, entity type 128, contain the only type of surface data which may be used directly. The control 
points and other data for these entities can be read directly from the IGES file. For other surface entity 
types such as planes, parametric splines, surfaces of revolution, tabulated cylinders, and ruled surfaces, the 
data must be converted into a NURBS or rational Bezier patch format. In all cases the parameter conversion 
results in a surface which is mathematically equivalent to the original surface. 

The second type of information used by RAGGS is information about curves. Curve entities may be in 
the form of parametric splines, NURBS, or point data. The curve data may be used subsequently for trim 
curves, surface generation, (surfaces of revolution, tabulated cylinders, ruled surfaces, etc.) or they may have 
other uses. Surface generation using point data involves curve fitting. This may result in a surface which is 
not exactly equivalent to the surface used to generate the points. 

Trimmed surface data is the third type of information which is retained and processed for later use. 
When combined with the surface and curve data, this data may be used to create the information necessary 
to define the active and inactive domains on each patch. This allows the trimmed surfaces to be displayed 
and to create grids which cross the trim boundaries. 

Ruled surface information is the fourth type of data which is stored. This data is combined with the 
curve data to create additional surfaces parameterized in the form of Bezier patches or NURBS. 

A file containing the patch information and the related boundary data read by RAGGS is referred to as 
a RASCAL file. 


NUMERICAL TOOLS 

Various numerical tools which are fundamental to our grid generations are described herein. 

Projecting A Point On A (Trimmed) Surface 

A point on a quilted surface may be identified by (p, s,f) with its position vector computed by i?(p,s,f). 
p is the patch number. For a given patch number p, s and t are the patch parametric variables and R is a 
Bezier patch or NURBS taking the form of Eq. (6). Let f denote a space point. Then, the square of the 
distance between R and r*is given by 

d= [fl(p,.s,t) - r] • [.R(p, a, f) - r] (16) 

The projection of f is defined to be the point R that minimizes d of Eq. (16). 


A local maximum or minimum distance is found by taking partial derivatives of d with respect to s and 
t as follows: 


R a * [r - R\ = 0 = /(p,M) 

R t ■ ff - ^1=0 = g(p,s,t) 
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(17) 

(18) 



Here, the subscripts s and t denote the derivatives with respect to s and t, respectively. Applying Newton’s 
iteration method, we obtain 

/,As + / t At=-/ ( 19 ) 

g a As + <7t At = -g ( 20 ) 


where As and At are the changes in s and t, respectively. Equations (19) and (20) are solved for (As, At) 
at fixed p. (s, t) are then updated as follows: 


( s Q + [w(|As| - c) + e] sign(As) if |As| > e 

\ s« + As if |As| < c 

f t^ + M| At | - e) + f] sipn(At) if |At| > <r 

\ t« + At if |At| < e 


( 21 ) 

( 22 ) 


Here, u; is the relaxation factor, c (typically, 0.02) is the cut-off parameter and the superscript Q denotes the 
Q-th iteration. 


During iteration, the patch number p is treated as a parameter and updated only if (s,t) is outside the 
active region of patch p. When the calculated (s,t) is outside the active region, the patch boundary data 
and the trim curve boundary data (if any) are used to find the adjacent patch to update p and (s, t) for the 
next iteration. In the case of patch boundary, Eq. (13) is used. In the case of trim curve, Eq. (15) is used 
to compute a b on the trim curve of the adjacent patch. The physical space point computed at <7 t on patch 
p is projected onto the trim curve on the adjacent patch using the computed a b as the initial guess. The 
corresponding (s,f) on the adjacent patch is determined by Eq. (10). Then the computed (s,t) is used to 

continue the iteration. 


Projecting A Point On A Curve And Distance Between Curves 

Consider two parametrically defined curves, f(t) and £(r), each taking the form of Eq. (1). The square 
of the distance between f(t) and f(r ) is given by 

d = [r(t) - f(r)] • \r{t) - f(r)] (23) 


At any point f(r) on the second curve we define the separation between the curves as the minimum distance 
d of Eq. (23). This is found by taking partial derivative with respect to t while fixing r as follows: 

n(t) ■ [r(0 - ((r)} = 0 = f(t, r) (24) 

Applying the Newton’s iteration method, we obtain 

f t At = -f (25) 


or 


A t = - 


n [r - f(r)] 


ft • r t + f tt ■ [f - f(r)] 


(26) 


where At is the change in t. 
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Along the curves there are a series of extrema where the curves are locally separated by a maximum or 
minimum distance. At these points, the partial derivative of Eq. (23) with respect to r vanishes at fixed t 
as follows: 

>r(*(r)) - 6(r)] • [m - f(r)] = 0 = f(t{r), r) (27) 

where t(r) is defined by the solution of Eq. (24) at fixed r. 

Applying the Newton’s iteration method, we obtain 

frAr = -f (28) 


or 


A r = 


~(n*r - Ct) • (r - f) 

{r t t TT + f u t\ - ( TT ) ■ (f- () + (f t t T - £ r ) • (f t t T - £ T ) 


Here, t T and t TT may by evaluated from Eq. (24). For example, t T is evaluated as follows: 


(29) 


t r 


ft ' (t 

rtt • (r - {) + f t ■ f t 


(30) 


So we can choose a r, iteratively solve for t and r from Eq. (26), and then get A r from Eq. (29). This 
process continues until Ar = 0. 


Projecting A Point On A Surface Curve 


Given a point in space f and a point on a parametrically defined curve, _ft(s(r), t(r)), on a surface R(s, t), 
the square of the distance between R and r is given by 


d = [j?[a(r), <(r)] - fj • [r[s(t), t(r)] - fj 


(31) 


A local maximum or minimum distance is found by taking partial derivative with respect to r as follows: 

R r • [5 - fj = 0 = f(r) (32) 


Applying Newton’s iteration method, we obtain 

fr At = -/ 


(33) 


or 

Ar- „ 

R t ' R r 4~ Rtt * (-ft — r) 

where 

Rj- — RgS-j- “I - RR't 

Rtt ~ Rs^tt 4" R'R , t t 4" 4" 4“ 4“ 


(34) 


(35) 

(36) 
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Curve And Surface Intersections 


Consider a curve C(u), and a surface R(s,t). At an intersection between C and R we have 

C(u)-R(s,t) = 0 = f(u,s,t) 


( 37 ) 


constituting 3 equations in 3 unknowns, u,s, and t. Following the Newton’s iteration method, we derive an 
iterative process as follows: 

-C u (u) Ait + R,(s, t)As + R t {s, t)At = C{u) - R(s, t ) (38) 


Isolating terms with cross and dot products, we have 3 equations, 

- R,(s , t) • [. R t (s , t) x C u (u)] Au = R a {s, t ) • {R t (s, t ) x [c(n) - R{s, t)] } 
C u (u) • [£ t (s, t ) x R a (s, *)] As = C u (u) ■ {R t (s, t ) x [(?(«) - R(s, t)] } 
C u (u) • [. R,(s , t) x R t (s, 0] At = C u (u) ■ { R,(s , t) x [c(u) - R(s , <)] } 


Solving these equations, we have the iterative process, 


Au = 


R a • [R t x /] ^ C u ■ [fit x /] ^ _ C» • [R a x /] 


* n 


C u • n 


, At = 


C u • n 


(39) 

(40) 

(41) 


(42) 


since 

Cu ■ [£, X it] = -Cu ■ [£f X R a ] = R a ■ [£ t x Cu] = Cu ■ n (43) 

where n = R a x Rt is the surface normal at (s,t), and f = C — R. For a quilted (trimmed) surface, patch 
number p will enter the iterative process and play the same role as that in section Projecting A Point On 
A (Trimmed) Surface.” 


Surface And Surface Intersections 

When surface intersections are not given by trim curves, an intersection must be computed. Two surfaces 
generally intersect along a curve or curves in space. Initial points on any intersection curves may be found 
by creating a mesh of curves, usually constant parameter curves, on one of the surfaces, and then finding the 
intersection of these curves with the other surface. The mesh must be of sufficient density to find at least 
one point on all curves of intersection, but not too dense and causing an inordinate number of computations. 
When a point on any intersection curve is found, further intersection points are easily created along this 
curve. This curve may then be described parametrically by fitting an interpolating function through the 
generated points. This procedure is followed for each curve of intersection. 

Consider two surfaces given by the equations 

£1 = P(s,t) (44) 

Ri = Q(u,v) (45) 

each taking the form of Eq. (6). At an intersection point, we have 3 equations in 4 unknowns s,t,u, and t>. 

P(s,t)~ Q(u,v) = 0 = f(s,t,u,v) (46) 
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(47) 


Following Newton’s iteration method, we have 

P 3 (s, t)As + P t (s, t)At - Q u (u, v)Au - Q v (u, v)Av - -f 
Since the system is overspecified, we need another equation for a unique solution. 


Holding one of the paramters constant reduces the problem to a curve and surface intersection. Holding 
a parameter constant is also useful for finding a point on a particular intersection curve, such as the patch 
boundary. We could also hold the parameter constant which maximizes the denominator of Eq. (42) for the 
curve and surface intersection. Another way is to approach perpendicularly to the intersection curve along 
one of the surfaces. The tangent vector to the intersection curve r, is given by, 

f = np X Hq (48) 

where Hp and Hq are the local normals at (s,t) and (u, v), respectively, 

n P - P s x Pf, Hq = Q u x Qv (49) 

The movement vector on surface P is given by P s As + P t At. The trace of the movement may be written in 
parameterized form as follows: 

C s As — P s As + P t At (50) 


with the direction of movement constrained to the normal of the intersection curve, i.e., 

f • (P s As + P t At) = 0 (51) 

Now, the problem has been reduced to the intersection of curve C(s) and surface Q(u,v). Following section 
” Curve And Surface Intersection” and from Eq. (42), we have 


As — — 


Qu ■ [q v x /] 


Au — — 


c, ■ [Qv x /] _ c 3 • [q u x f\ 


Av = 


* riQ 


C 3 * UQ * n Q 

Here, C s may be calculated from Eqs. (50)-(51) as follows: 

(f-P t )C s As = (f • P t )(P a As + P t At) 

= [(T ■ P t )P, - (T ■ P a )P t ] As 

(f • P t )P 3 - (f ■ P a )P t = [(n P x uq) x n P ]A« 

Since 

(f • P t )C a ■ Uq = [(np X Uq) X np] • Uq — (np X Uq) ■ (np X Uq) 


(52) 


(53) 

(54) 

(55) 


we cannot have a nonzero denominator unless rip x uq = 0, which means the surfaces are locally parallel 
and there is no intersection anyway. Substituting in the above we get the iteration process, 


As = 

[(np x n Q ) ■ P t ] [<9„ x /] • Q u 
(np x Hq ) * (np x wq) 

(56) 

At = 

[(np X Uq) ■ P„] [Qv x /j • Qu 

(57) 

(rip x uq) • (np X Hq) 

Au = 

[(np x ng) x np] • [<?„ x /] 
(np x ng) • (np x ng) 

(58) 

Av = 

[(np X ng) X np] • [<5„ X /] 

(59) 

(np x Hq) • (np x Hq) 
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As previously defined, ftp and nq are in the direction of the local normals. 

np = P,x P t , nq = Q u x Q v , /= P(s,t) - Q(u,v ) 


(60) 


When the convergence is satisfied, the next point may be found by marching along the surfaces in the 
direction of the tangent vector f = n P x nq, and iterating again. Along the tangent vector, 


A R = P 3 As + P t At = Act 

In 


(61) 


or 


A R = Q u Ati + Q v Av = — tAct 

In 


(62) 


where a is the arc length. Then As and At along the intersection curve may be calculated by the vector 
cross product of Eq. (61) and P 3 and P t , respectively, 


—(Pa X P t )As 

As 

At 


P t x [(A x P t ) x n Q ] = (P, x P t )(Pt ■ ”Q)^ 

- - Act 

~(Pt ■ n Q)-j^j- 

/ f? - . Act 

(P,-nq)-^ 


(63) 

(64) 

(65) 


Similarly, Au and Av are obtained by the vector cross product of Eq. (62) and Q u and Q v , respectively, 


A u = (Q v ■ np)—— 

l r l 

Au = -(Q u • ” p )-|fj- 


( 66 ) 

(67) 


Here, A t is computed as follows: 

At = t s As + f t At + f u Au + f v Av 

= [~T,(P t ■ nq) + f t (P 3 ■ nq) + t u (P v • np) - t v (P u ■ np)] -pj- 


( 68 ) 


Surface And Plane Intersections 

For the user’s convenience, a space plane is specified by a set of three non-colinear points or by a reference 
point and a normal to the plane. Intersection between a quilted surface and the user-specified plane is done 
exactly in the same manner as that described in section ” Surface And Surface Intersections” once the plane 
is represented by a quilt containing a single patch. 

SURFACE GRID GENERATION 

Both structured and unstructured grid generators used in RAGGS are briefly described in this section. 
References will be cited for readers who are interested in details of the numerical algorithms actually imple- 
mented for the surface grid generations. 
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Structured Surface Grid Generation 


TFI Method. — The transfinite interpolation (TFI) of Ref. 2 is coupled with the space point projection, 
Eqs. (16)-(22), to determine the interior grids of a 4-sided block. A subsection approach is employed. The 
approach allows a 4-sided block grid be calculated subsection by subsection. Each subsection is defined by 
its two opposite diagonal points and can easily be specified interactively by point- and- click. A subsection 
grid is computed in two steps. First, the 4-sided subsection grid is calculated by TFI based on the sectional 
edge grids. Second, Eqs. (19)-(20) are used to project the TFI grids onto the quilted surface. At each point, 
a set of initial guesses (p, s,<) is formed by gathering all the available values at immediately adjacent grid 
points. Additional sets of (p, s,t) may also be specified by user. The solution is taken to be the one that has 
the minimum distance among all the calculated solutions. 

When a grid is generated on a single patch, TFI may be optionally applied to the surface parametric 
variables. In this case, the grid position vectors are evaluated directly from the patch considered using the 
TFI interpolated values of parametric variables. 

PDE Method. — The elliptic grid generation system on arbitrary curved space surface can be written, in 
tensor form 3 , as follows: 

Y Y + Y 99 nn P n rt" - gHn = 0 (69) 

m = l n=l n=l 

Here r is the grid position vector, (* are the curvilinear coordinates, P n are the grid control functions, g il 
is the contravariant metric tensor, and g = det(gij) with gij the covariant metric tensor. H is the twice 
of the surface mean curvature and n is the unit surface normal. Equation (69) is solved numerically using 
the relaxation methods described in Ref. 4 once the differential operators are replaced by central- difference 
operators. Equations (19)- (20) are used to snap the solution grid points onto the quilted surface at the end 
of each iteration. 

The grid control functions P n are evaluated at the 4-sided block edges and the boundary values are inter- 
polated into the field. Both Thomas/Thomas-Middlecoff’s method 5,6 and Sorenson’s method 7,8 have been 
implemented for the construction of orthogonal boundary grids. However, Sorenson’s method is implemented 
based on Thompson’s method of Ref. 2 with the following two modifications. One, the central-difference 
operator is used to approximate the second-order derivatives normal to the boundary instead of the one-sided 
difference operator. Second, the grid control functions are updated immediately after the contributions from 
a given orthogonal boundary section have been evaluated and the control function increments are evalu- 
ated against the most current values. The first modification provides a consistent numerical method in the 
sense that both the grid position vector and control functions satisfy the same governing equation as the 
solution converges, and a grid initially orthogoanl will remain and be returned as the solution. The second 
modification has provided a more stable iterative process. 

Similarity Method. — An existing grid may be fine-tuned to have a similar or identical distribution in a 
specified base variable, the arc length or one of the three Cartesian coordinates, along a specified coordinate 
direction. To accomplish this, the edge distributions of the specified base variable are interpolated into the 
field. Then, the grid points are re-evaluated line by line in the specified coordinate direction and projected 
onto the quilted surface. This can provide a simple means of achieving a cylindrical-coordinate-like grid. 
Also, it can provide panels on an aircraft configuration with constant spanwise or axial sections as required 
by some panel codes. 

Marching Method. — A grid may be constructed on a quilted surface from a given line grid (open or 
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closed) by marching in the direction normal to the line segment with a specified step size The marching 
direction is the direction of the vector formed as the cross product (in the order specified) of the surface 
normal and the tangent vector to the line segment. 

Smoothing.— An existing grid may be smoothed by replacing each of the interior points by the averaged 
value of its adjacent grid points. The number of passes of smoothing is specified by the user. 

Unstructured Surface Grid Generation 

Details of our unstructured grid generation algorithm may be found in Ref. 9. Briefly, the technique 
is based closely on the Lohner-Parikh’s 10 data structures and the 2-D advancing front triangulation sc erne 
described in Ref. 11. However, the implementation is made in the 3-D physical space coupled with 
space point projection, Eqs. (16)-(22). Interactively, a grid generation starts with a rectangular octree box 
specified just big enough to cover the entire region to be triangulated. The advancing front is constructed 
automatically as the user is pointing and clicking line segments on the screen. During triangulation, every new 
nodal point is introduced to form a nearly equilateral triangle and satisfy the mesh consistency requiremen 
with a base selected from the front. Equations (19)-(20) are used to iteratively locate the new nodal point 

on the quilted surface. 

Since the unstructured grid generation is based on the same geometry database as that of the structured 
grid generation, geometry manipulators and domain blocking tools developed for structured grid generation 
can be used directly. With this new capability, RAGGS provides a unique environment where both structured 
and unstructured grids can be generated and manipulated on a common geometry database. 

SAMPLE CASES 

Sample cases are presented herein to demonstrate the capabilities of RAGGS in handling trimmed surfaces 
and generating grids on an unstructred quilt of patches and/or trimmed surfaces. 

Figure 4 demonstrates the capability of RAGGS to display trimmed surfaces. Shown in the figure is a 
portion of an impeller. The geometry is defined by a series of trimmed surfaces. Figures 4-(a) and (b) show 
the shaded displays of the untrimmed and trimmed patches, respectively One can hardly picture what 
geometry would look like from the untrimmed surface display in Fig. 4-(a). But it is clearly shown g_ 
4-fb) Figures 4-(c) and (d) show how the two patches corresponding to the front face and the right side of 
the geometry are discretized for display purposes. In RAGGS a structured grid is used to discretize the patch 
in the parametric space with the grid size adaptive to the local curvature in each direction of the parametric 
variable as shown in Fig. 4-(c). For a trimmed patch, the active region is automatically triangulated in 
parametric space, using the advancing front technique of ReL 9, with the trim curves as the ini i tot ^ 
Figure 4-(d) shows the triangulation of the active regions of the two patches shown in Fig. 4-(c) It may be 
noted that in RAGGS the trimmed and untrimmed displays can easily be toggled on and off. Also noted is 
that no additional patch definition is given to the trimmed patch and the definition of the original untnmmed 

patch is used for the geometry. 

Figure 5 demonstrates the capability ot RAGGS for generating structured and unstructured grids on 
unstructured quilts of patches. Figure 5-(a) shows a B-1B forebody with fins Four quilts are shown here: 
the upper and lower body quilts and the upper and lower fin quilts. No trim curves were provided so a 
surface intersection was computed between the body and fin qnilts. The irregular arrangement of patches is 
clearly indicated by the wireframe display in Fig. 5 (a). The calculated structured and unstructured grids 

are shown in Figs. 5-(b) and (c), respectively. 
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Figure 6 demonstrates the capability of RAGGS for gridding on trimmed surfaces. Figure 6-(a) shows a 
plane and a bump (part of semi-spherical surface represented by NURBS). The four edge grids are specified 
entirely on the flat plane, and the interior grids are generated while considering the trimmed surfaces. As 
shown in Fig. 6-(b), the interior grid points have been properly positioned on the trimmed patches. It is 
noted that the unstrimmed surfaces are the only surfaces used to define the geometry. 

CONCLUSIONS 

Capability of generating structured and unstructured grids directly on CAD generated geometry with/whitout 
trimmed surfaces has been described and demonstrated. Both geometry integrity and patch structure are 
retained from IGES file to CFD geometry database. 
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Figure l.-The program structure of RAGGS. 



Figure 2-Typical boundary connectivity between two adjacent patches and 

determination of new patch and parametric variables during iteration. 
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(b) Shaded display of trimmed surfaces W) Triangulation ot trimmed patches 


Figure 4. -Capability of displaying trimmed and untrimmed surfaces. 
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(b) Structured surface grid 


Figure 5.-Capability of structured and unstructured grid generations on unstructured quilt of patches. 
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(a) Patch-defined plane and bump 



(b) A structured grid on the trimmed surface 


Figure 6. -Capability of gridding on trimmed surfaces. 
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